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Abstract 

Nuclear propulsion has been identified as 
one of the key technologies needed for human 
exploration of the Moon and Mars. The 
Nuclear Thermal Rocket (NTR) uses a nuclear 
reactor to heat hydrogen to a high temperature 
followed by expansion through a conventional 
convergent-divergent nozzle. A parametric 
study of NTR nozzles was performed using the 
Rocket Engine Design Expert System 
(REDES) at the NASA Lewis Research Center. 
REDES used the JANNAF standard rigorous 
methodology to determine nozzle performance 
over a range of chamber temperatures, chamber 
pressures, thrust levels, and different nozzle 
configurations. A design condition was set by 
fixing the propulsion system exit radius at five 
meters and throat radius was varied to achieve 
a target thrust level. An adiabatic wall was 
assumed for the nozzle, and its length was 
assumed to be 80% of a 15° cone. 

The results of this study conclude that 
although the performance of the NTR, based on 
infinite reaction rates, looks promising at low 
chamber pressures, finite rate chemical reac- 
tions will cause the actual performance to be 
considerably lower. Parameters which have a 
major influence on the delivered specific 
impulse value include the chamber temperature 
and the chamber pressures in the high thrust 
domain. Other parameters, such as two-dimen- 
sional and boundary layer effects, kinetic rates 
and number of nozzles, affect the deliverable 
performance of an NTR nozzle to a lesser 
degree. For a single nozzle, maximum perfor- 
mance of 930 seconds and 1030 seconds occur 
at chamber temperatures of 2700 and 3100 K, 
respectively. 
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Introduction 

Nuclear propulsion has been identified as 
one of the key technologies needed for human 
exploration of the Moon and Mars. 1 Successful 
application of any mission using an NTR 
depends on high delivered performance, 
because the weight of a nuclear reactor must be 
compensated for by reduced propellant require- 
ments. To date, many different reactor design 
and operating conditions have been identified 
as a means of satisfying various mission 
scenarios. The different reactor operating 
conditions include both pressure and thrust 
range variations of nearly two orders of magni- 
tude, and reactor temperatures operating up to 
the melting point of the fuel. 2,3 With such a 
variety of perceived operating conditions, a 
priori knowledge of anticipated nozzle perfor- 
mance at a certain reactor operating condition 
is important for design and evaluation 
purposes. Likewise, to design the optimum 
NTR system, the sensitivity of nozzle perfor- 
mance to reactor operating conditions is 
crucial. 

The purpose of this paper is to investigate 
the effect of the various operating conditions on 
the NTR nozzle performance in a rigorous 
manner, and to identify the major parameters 
which affect the deliverable performance of an 
NTR nozzle. In order to determine their affect 
on the delivered nozzle performance of an 
NTR, a variety of operating parameters were 
investigated over a wide range of nozzle 
configurations. Chamber pressure (p c ), 
chamber temperature (T c ), propulsion system 
thrust level (F), and number of nozzles were 
investigated with four different combinations 
of nozzle performance codes to determine their 


Table 1. Analysis parameters and values. 


Parameter: Values 

Chamber Pressure : 

10, 30, 100, 300, 1000, 3000 psia 

Chamber Temperature : 2700, 2900, 3100 K 

Propulsion System Thrust Level : 

10,30, 100, 300 Klbf 

Number of Nozzles : 1, 3, 5, 7 


throat radius and area ratio of the nozzle was 
allowed to vary to meet the target thrust level 
constraint. 

Analysis and Results 
Results of the studies which describe the 
Figure 1. Mole fraction of H 2 for different 
chamber pressures at T c =3100 K. 

0.50 


effect on specific impulse (1^). Ranges over 
which the parameters were varied are given 
in Table 1. In each case, rocket nozzle 
performance calculations were made to 
determine the chemical equilibrium perfor- 
mance, and the detrimental effects of finite- 
rate chemistry, two-dimensional flow, and 
boundary layer growth. A matrix was 
investigated using a computer code which 
conforms to the Joint Army -Navy -NASA- 
Air Force (JANNAF) standardized method- 
ology for determining chemical liquid 
rocket nozzle performance. 

A nozzle cluster exit radius of 5 meters 
was used throughout the analysis. This 
radius corresponds to the payload bay size 
of the advanced launch vehicles formerly 
and currently under consideration for 
development (e.g., Shuttle-C and ALS). 4 
Using this criteria, the entire nozzle 
assembly could fit into the cargo bay of 
these vehicles. As will be shown, the use of 
this maximum practical nozzle cluster size 
as a design criteria corresponds to a nearly 
maximum performance condition. 

In order to generate design study results 
with fixed exit area and thrust level as a test 
matrix parameter, a need existed to over- 
come the limitation of available codes 
which calculate nozzle thrust based on 
given input parameters. This was accom- 
plished through the use of the Rocket 
Engine Design Expert System currently 
under development at the NASA Lewis 
Research Center. 5 At each data point, the 
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recombination of monatomic to diatomic 
hydrogen are presented and discussed first, 
since this phenomenon is the prime contributor 
to the NTR nozzle’s high performance. Next, 
the effect of varying the operating parameters 
on nozzle geometry and performance are given. 
How performance is affected by utilizing 
multiple nozzles in a cluster is discussed next. 
Finally, the uncertainty of chemical recombina- 
tion rates on nozzle performance is given. 

Hvdrogen Recombination 

The energy release which accompanies the 
recombination of monatomic hydrogen (H) to 
diatomic hydrogen (H 2 ) provides the high 
performance for which the low pressure NTR 
concept is known. Understanding the recombi- 
nation process is key to understanding how the 
NTR performs. 

Mole fractions of H 2 , computed using one- 
dimensional infinite and one-dimensional finite 
chemical reaction rates, were plotted as a func- 
tion of subsonic and supersonic area ratios for 
different thrust levels and p c values (Figure 1). 
In all cases, the H 2 mole fraction starts at a 
minimum value in the chamber. The infinite 
and finite reaction rate results in this part of the 
chamber are identical because the residence 
time of the propellant is much greater than the 
time required for the reactions to reach equilib- 
rium. thereby allowing the recombination reac- 
tions to reach equilibrium. As the flow 
expands, moving in Figure 1 from left to right 
along the curve toward the throat (where the 
area ratio equals one), and continuing along the 
supersonic portion of the curve, the rise in the 
H ; mole fraction indicates the recombination of 
H into H 2 . In all cases, it is seen that the larger 
thrust nozzles are closest to being in chemical 
equilibrium. This occurs because the larger 
thrust nozzles also correspond to larger throat 
dimensions, and a more gradual fluid accelera- 
tion rate. The more gradual acceleration rates 
allow appreciably more H recombination to 
occur, pushing performance toward the infinite 
reaction rate limit. 


Superimposed on the recombination plots, 
in dotted lines, is the first derivative of the mole 
fraction curves. Where recombination occurs, 
the derivative is denoted by a non-zero value. 
Frozen conditions occur when the derivative is 
zero in the supersonic portion of the nozzle. As 
can be seen in Figures la, lb, and lc, the 
majority of H 2 recombination and its associated 
energy release occurs in the near throat region. 

Nozzle Geometry 

A capability of the Rocket Engine Design 
Expert System (REDES), called the Thrust 
Dependent Engine Sizing (TDES) function 6 , 
iterated toward a desired thrust level while 
employing the rigorous JANNAF solution 
procedure 7 to calculate nozzle 1^. A liquid 
rocket nozzle performance evaluation code, 
Two Dimensional Kinetics (TDK) 8 , was used 
for the entirety of this analysis. 

Table 1 gives the parameters included in the 
study and their values. Geometric parameters 
which define the contraction and near-throat 
region of the convergent-divergent nozzle are 
given in Table 2. Other assumptions which 
were used to characterize the nozzle shape 
included: 

• the nozzle contour was a parabola 
fitted to connect a tangent point 
on the downstream throat circle 
and a specific exit coordinate, 

• the exit radius of the entire pro- 
pulsion system package was set to 
be 5 meters, 

• nozzle length was determined by 


Table 2. Baseline geometric parameters. 


Chamber Contraction Ratio = 

5.00 

Nondimensional Throat Upstream Radius of 

Curvature = 

2.0623 

Nondimensional Throat Downstream Radius of 

Curvature = 

0.1900 

Throat Upstream Tangent Angle = 

26.25” 

Throat Downstream Tangent Angle = 

32.00* 

Nozzle Length = 80% of a 15” Cone 
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Figure 2. Nozzle area ratio variation as a function of chamber 
pressure, for T c =3100 K, one nozzle. 



the TDES function of REDES to 
be 80% of a 15° cone nozzle with 
the same throat radius, down- 
stream radius of curvature, and 
area ratio, and is calculated from 
the following equation, 9 

1) +RWTD ( seca- 1) 


where L n is nozzle length, r, is the 
throat radius, £ is the area ratio, 
RWTD is the downstream throat 
radius of curvature, and a is the 
cone half angle (15° in this case), 
and 

• there was no spacing between 
nozzles in the 3,5, and 7 nozzle 
cases. 

Boundary layer assumptions included: 

• an adiabatic wall, 

• equilibrium chemistry in the 
boundary layer, 

• the standard JANNAF boundary' 
layer loss methodology was used, 
and 


• the boundary layer would become 
turbulent after the flow Reynold’s 
number based on momentum 
thickness achieved a preset value 
of 360. 

Because the TDES function within REDES 
kept the nozzle exit radius constant and varied 
the throat radius to arrive at a desired thrust 
level, the area ratio of the nozzle varied as a 
function of chamber pressure and thrust level. 
Figure 2 shows a typical plot of nozzle area ra- 
tio versus chamber pressure for all thrust levels 
for a chamber temperature of 3100K and the 
single nozzle case. All plots of nozzle area ra- 
tio as a function of chamber pressure are given 
in Appendix. 

As can be seen in Figure 2, nozzle area ra- 
tios of 60,000:1 result in some instances. This 
is due to the combination of low thrust, high 
chamber pressure, and the fixed exit area con- 
straint. In practice, relaxation of the fixed area 
constraint would likely occur because of noz- 
zle weight and performance considerations. 
Small area ratios, less than 10:1, occur at high 
thrust, low chamber pressure operating condi- 
tions. 

Figure 3 shows the nozzle length trend for 
the same chamber temperature and number of 
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nozzles as Figure 2. All nozzle length plots are 
given in Appendix. Nozzle lengths at low p c 
levels are very sensitive to thrust level. Nozzle 
lengths for higher thrust levels are shorter than 
those for lower thrust levels, because of the 
lower area ratios for those nozzles. The differ- 
ence, however, between the 10Klb f and the 
300Klb f nozzle lengths decreases drastically 
and the nozzle lengths converge toward a sin- 
gle value as chamber pressure increases. This 
is the result of the area ratio increase with 
chamber pressure, due to the fixed exit area de- 
sign condition. As area ratio increases, the 
nozzle length equation approaches a value of 
approximately 3r e (shown as the horizontal 
dashed line in Figure 3). For the single nozzle 
case, values converged toward a length of ap- 
proximately 590 inches at high chamber pres- 
sures. In the case of three nozzles, the 
asymptotic length value is roughly 270 inches, 
for five nozzles, lengths approach 210 inches, 
and in the seven nozzle case, lengths converge 
toward a value of 200 inches. 

Nozzle Performance 

Equilibrium performance of a NTR nozzle 
can be predicted to exceed 1200 seconds as 
seen in the case presented in Figure 4. 


However, equilibrium performance of a NTR 
nozzle must be decremented by real effects and 
performance losses which can be predicted 
using the JANNAF methodology, standardized 
nozzle performance code, TDK. 8 TDK calcu- 
lates performance decrements due to chemical 
kinetics, two-dimensional losses, and boundary 
layer growth. Other loss mechanisms, 
including mixing, vaporization and non- 
boundary layer heat losses, were assumed to be 
negligible, or not applicable to a nuclear 
thermal rocket nozzle simulation. After 
subtracting these losses from the equilibrium 
specific impulse (1^) values, the resulting 
performance was the delivered 1^. 

Results including all loss mechanisms were 
plotted to show NTR nozzle performance for 
the entire range of chamber pressures, chamber 
temperatures, and thrust levels. Figure 5 is a 
representative plot from the total set which 
were created over the range of parametric vari- 
ables (Table 1). All plots are given in the 
Appendix. 

To identify the magnitude of each loss 
mechanisms, nozzle performance was 
computed using each module of the TDK 
program in a sequential fashion. Isp values 
based on one-dimensional infinite chemical 


Figure 3. Variation of nozzle length as a function of chamber 


pressure. 
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Figure 4. Theoretical equilibrium chemistry performance of a 

nuclear thermal rocket nozzle with exit radius of 5 meters. 

Chamber temperature = 3100 K, thrust level = 10,000 Ibf, 


single nozzle case. 



reaction rates were calculated with the One- 
Dimensional Equilibrium (ODE) subprogram. 
Next, the one-dimensional finite chemical reac- 
tion rate performance predictions were 
computed with the One-Dimensional Kinetics 
(ODK) subprogram, followed by an inviscid 
Two-Dimensional Kinetics (TDK) subprogram 
calculations. Finally, the Boundary Layer 
Method (BLM) subprogram, with a final TDK 
run, was added to this series. The second pass 
of TDK calculated the nozzle performance in 
the boundary 7 layer-displaced inviscid core. 
Therefore, I sp values were calculated by 
running the series of ODE, ODK, TDK, and 
TDK -B LM -TDK codes. When plotted 
together, the gaps between the curves labelled 
ODE and ODK, and between ODK and TDK 
show how chemical kinetic and two-dimen- 
sional effects vary as a function of chamber 
pressure, respectively. The gap between the 
TDK and the TDK-BLM-TDK curves show the 
variation of boundary layer losses with 
chamber pressure. 

The two-dimensional chemical kinetic and 
boundary’ layer losses taken together are accu- 
rate. However, due to coupling of the two- 
dimensional and boundary layer calculations, 


the isolated two-dimensional and boundary 
layer losses can only be considered approxi- 
mate. Isolated boundary layer losses must be 
interpreted carefully since the TDK curves are 
plotted from data computed during the first 
TDK computational pass. The subsequent TDK 
pass, using a boundary layer displaced geom- 
etry, increases the inviscid core performance. 
To legitimately show the boundary layer drag 
losses as a function of chamber pressure, 
another curve would have to be plotted, based 
on the k. p calculated during the second pass of 
TDK in the TDK-BLM-TDK computational 
progression. 

Infinite reaction rate I sp (ODE) specific 
impulse start very high at low chamber pres- 
sures and decreases as p c increases. This is due 
to the increased levels of dissociated hydrogen 
in the chamber at lower chamber pressures, 
which contributes directly to the overall perfor- 
mance through the increased energy released 
when recombination occurs, and the perfor- 
mance also increases due to the decreased 
molecular weight of the propellant. As chamber 
pressure increases, the recombination of H into 
H ; decreases since the initial amount of disso- 
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Figure 5. Specific impulse breakdown for 3100 K chamber 

temperature and 10KJb f thrust level, 1 nozzle with exit 
radius of 5 meters. 
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dated hydrogen in the chamber is diminished at 
higher chamber pressures. Only in one case 
(T c =3100 K, and F=300,000 lb f , one nozzle, 
Figure A- 1 6) does the equilibrium performance 
fail to follow the trend described above. This is 
due to the very large throat required to produce 
the desired thrust at the given chamber pres- 
sures. A large throat results in a lower area ratio 
(< 1 0: 1 ), which reduces the performance of the 


nozzle. In all other cases, nozzle area ratio was 
a secondary consideration compared to 
hydrogen recombination effects. 

As can be seen in Figure 5, the ODE perfor- 
mance results are greater than the ODK values 
at all chamber pressures, although the differ- 
ence is greatest at lower chamber pressures. 
Because of relatively slow chemical recombi- 
nation at low chamber pressures, ODK perfor- 


Figure 6. Decrements to I due to one-dimensional kinetic 

effects as a percentage of ODE performance for all 
operating conditions. 
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Figure 7. Kinetic performance efficiency (I^ODK) divided by 
I sp (ODE)) depicting kinetic losses for all operating 
conditions. 



mance, I sp (ODK), is appreciably less than 
I sp (ODE). Figure 6 depicts the kinetic losses as 
a percentage of equilibrium 1^. As chamber 
pressure increases, the ODE and ODK values 
converge in Figure 5, because as chamber pres- 
sure increases, the amount of H present is 
minimal and because recombination rates 
increase significantly as a function of pressure. 


This trend is reflected in Figure 6 by the kinetic 
loss curves moving closer to the horizontal 
axis, which represents very small kinetic 
losses. 

Kinetic losses of 15% occur at low chamber 
pressure and thrust level, the precise operating 
conditions where the equilibrium specific 
impulse was predicted by ODE to be in excess 
of 1200 sec (refer to Figures 4 or 5). This is the 


Figure 8. Performance efficiency (IjpCTDK-BLM-TDK) divided 
by Isp(ODE)) depicting kinetic, two-dimensional, and 
boundary layer losses for all operating conditions. 
Single nozzle case. 
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major reason why such high nozzle perfor- 
mance cannot be attained at these operating 
conditions. 

Kinetic performance efficiency, defined as 
the IjpfODK'' divided by the I^ODE), is shown 
for the entire matrix of p c , T c , and F in Figure 7. 
As the extent of dissociated hydrogen present 
in the chamber increases monatomically with 
increasing chamber temperature, the overall 
trend of decreased performance efficiency with 
increasing chamber temperature would be 
expected, as shown in Figure 7. 

Two-dimensional and boundary layer 
losses, as computed in the TDK-BLM-TDK 
sequence of analyses, account for the perfor- 
mance decrement between the curves labelled 
ODK and TDK-BLM-TDK in Figure 5. 
Combined, the two-dimensional and boundary 
layer losses subtract between 1.5 and 3.5% 
from the I sp (ODE). Compared to a possible 
15% loss due to kinetic effects, these losses 
seem minor. However, at high chamber pres- 
sures where kinetic effects are minimal, the 2% 
two-dimensional/boundary layer losses are the 
major contributor to performance degradation. 
Figure 8 shows the delivered performance effi- 
ciency, I sp (TDK-BLM-TDK) divided by 


I^p(ODE), which includes the boundary layer 
and two-dimensional losses along with the 
kinetic losses. This figure indicates that at low 
pressure, low thrust level, and for an exit radius 
of 5 meters, a performance efficiency of 0.82 
can be expected. This lowers the predicted 
Isp(ODE) of 1230 seconds to a delivered 
Isp (TDK-BLM-TDK) of 1010 seconds. 

The actual value of delivered 1^ at each 
point provides the most important glimpse as to 
which parameters most affect the delivered 
nozzle performance of an NTR. Delivered I sp as 
a function of chamber pressure, for all chamber 
temperatures and thrust levels, is shown in 
Figure 9. From this figure, it can be seen that 
maximum deliverable 1^ values range between 
930 seconds for a T c of 2700 K, to 1030 for a T c 
of 3100 K. 

Also from Figure 9, it can be seen that for 
every 200 K increase in T c , a 40 to 50 second 
increase of delivered 1^ results. At low 
chamber pressures, the thrust level can greatly 
affect the nozzle performance. Figure 9 also 
shows that, to get a maximum 1^ value at low 
chamber pressures, a low thrust level (corre- 
sponding to high area ratios) would be desir- 
able. For low chamber pressures, changing the 


Figure 9. Delivered nozzle performance, including kinetic, two- 
dimensional, and boundary layer losses of an NTR. 
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Figure 10. Cluster configurations of 3, 5, and 7 nozzles used in this study. 



r/R = 0.46 

Area = 64% of 1 Nozzle 
Three Nozzles 



r/R = 0.36 

Area = 65% of 1 Nozzle 
Five Nozzles 



r/R = 0.33 

Area = 78% of 1 Nozzle 
Seven Nozzles 


thrust level from 10,000 to 300,000 lb f 
decreased performance by approximately 80 
seconds. Higher performance at low thrust 
occurs because this condition corresponds to a 
high area ratio nozzle, which more than 
compensates for the reduced kinetic efficiency 
of a low thrust nozzle. However, at moderate to 
high chamber pressures (p c >100 psia), the 
selection of thrust level has little effect on 
performance. This indicates that I sp at high p c is 
insensitive to area ratio variations. 

Multiple Nozzles 

An investigation into the effect of multiple 
nozzles on the performance and geometric 
parameters was conducted. Clusters of 1, 3. 5, 
and 7 nozzles which fit into the same five meter 
exit radius were investigated. Figure 10 shows 
the arrangement of each multiple nozzle 
configuration. Ratios of the single nozzle 
radius to the multiple nozzle radii, as well as 
the percentage of area coverage by the multiple 
nozzle configurations, are given. In the study, 
each configuration had an equal propulsion 
system thrust level. 

For the most part, multiple nozzle clusters 
consisting of 3, 5, or 7 nozzles performed 
between 1 to 2% (10 to 20 seconds of I sp ) below 
their single nozzle counterparts. In only a few 


cases (low chamber pressure and high thrust 
level for all chamber temperatures) did the 
performance penalty of going to multiple 
nozzles exceed 2% and reach approximately 
4%. Because of this, the resulting performance 
curves for multiple nozzles closely resembled 
those shown in Figure 9. Delivered perfor- 
mance for T c =3 100 K and F=300Klb f is plotted 
for 1, 3, 5, and 7 nozzles in Figure 11. The rela- 
tive insensitivity of nozzle performance to 
multiple nozzle configurations was expected 
because: 

• boundary layer losses were shown 
to be small, 

• two-dimensional losses are insen- 
sitive to chamber pressure and 
thrust levels, and 

• little hydrogen recombination 
occurs beyond the throat. 

All multiple nozzle curves are given in Appen- 
dix, Figures A-19 through A-72. 

A benefit of the use of multiple nozzles is 
the geometric compactness of multiple nozzles 
as compared to single nozzle geometries. 
Figure 1 2 shows a comparison of nozzle 
lengths for 1, 3, 5, and 7 nozzles. As the 
number of nozzles (engines) increases to attain 
a fixed thrust level, the amount of thrust per 
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Figure 11. Effect of multiple nozzles on performance. Chamber 
temperature = 3100 K, 300,000 lb f thrust case. 



nozzle decreases, but so does the exit radius per 
nozzle. The throat radius computed to attain the 
target thrust level resulted in area ratio trends 
which follow the percentage of a single nozzle 
area trends (given in Figure 10). A 65% reduc- 
tion in nozzle length can be achieved by 
utilizing a five nozzle cluster instead of a single 
nozzle to produce the same thrust level. Gains 
attained by increasing the number of nozzles 
from 5 to 7 are less than 5%. The results of 
Figure 1 2 combined with the finding that 
overall delivered performance remained 
roughly constant with the use of multiple 
nozzles indicates that substantial reductions in 
the overall nozzle package can be realized. 
Simplified analyses of nozzle weights (based 


only on surface area) indicate that mass savings 
of approximately 25% can be made for the 
cluster of five nozzles. 

Therefore, although multiple nozzle config- 
urations may decrease slightly the nozzle 
performance, it is clear that shorter nozzles can 
be designed by using a cluster of multiple 
nozzles which, in turn, permits packaging of 
the propulsion system into a smaller volume. 

Chemical Reaction Rates 

The chemical reaction rates used in this 
study (given in Table 3) were based on sugges- 
tions made by the National Aerospace Plane 
Rate Constant Committee of the NASP High- 
Speed Propulsion Technology Team 10 . 


Figure 12. The effect of multiple nozzles on overall nozzle length, expressed as a percentage 



% Single Nozzle Length 
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Table 3. Baseline reactions and rates. 
Reaction: Reaction Rate Equation 

H + H + H, -» H 2 + H 2 with k = 

H + H + H— »H 2 + H with k = L5xlO^T^ 

A sensitivity study of the effect of chemi- 
cal reaction rate on specific impulse was per- 
formed. Chemical reaction rates are calculated 
in the TDK program by the Ahrennius equa- 
tion, 

1000b 

k = AT n e RT 

where k is the reaction rate, A, b, and n are rate 
constants, T is the reaction temperature, and R 
is the specific gas constant. To determine sensi- 
tivity of the results to the reaction rate constant 
values, a series of computer runs were 
performed for the 100,000 lb f thrust level, the 
2900 K chamber temperature, and the single 
nozzle case. Reactions which were an order of 
magnitude slower and faster than the baseline 
reaction rates were input and the performance 
was calculated using the ODE-ODK program 
sequence to determine the effect of reaction 
rate variation on kinetic losses. The uncertainty 
in H recombination rates is currently thought to 


be between a factor of 2 and 4. 107 By varying 
the baseline reaction rates by an extreme factor 
of uncertainty, 10, and recomputing the NTR 
nozzle performance at one set of operating 
conditions, the change in specific impulse for 
an order of magnitude change in reaction rate 
was determined. Results are shown in Figure 
13. 

For chemical reactions which were a factor 
of 10 slower than the NASP suggested rates, 
Isp values were approximately 10 seconds (1%) 
lower for chamber pressures under 300 psia. 
Above 300 psia, the effect of slowing down 
the reaction rates seemed to diminish. 

For chemical reactions which were a factor 
of 10 faster than the suggested rates, the trend 
of 1^ was approximately the same as for the 
slower reactions. An increase of approximately 
20 seconds (2%) was seen in the cases where 
chamber pressure is less than 300 psia. In the 
300 psia case and above, however, a smaller 
increase in performance was predicted. This is 
due to nearly negligible kinetic losses at high 
chamber pressures. Overall, the influence of 
chemical reaction rates at all chamber pres- 
sures is secondary, causing a maximum varia- 
tion of approximately 2%. 


Figure 13. Effect of recombination rate variation on NTR nozzle 
performance. Chamber temperature is T C =2900K, 
F=100Klbf. single nozzle case. 
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Summary of Results 

Nozzle performance of a Nuclear Thermal 
Rocket was calculated over a range of chamber 
pressures, chamber temperatures, and thrust 
levels. NTR performance delivered by clusters 
of multiple nozzles were also considered. 

Recombination of monatomic hydrogen to 
diatomic hydrogen was investigated. Results 
show that a majority of the recombination 
process occurs in the near-throat region. 

Nozzle lengths at low p c levels are very 
sensitive to thrust level. Nozzle lengths for 
higher thrust levels are shorter than those for 
lower thrust levels, because of the lower area 
ratios for the high thrust nozzles. As p c 
increases, nozzle lengths for all thrust levels in 
different multiple nozzle configurations 
asymptotically approach lengths of 3r e (1 
nozzle), 1.4r e (3 nozzles), l.lr e (5 nozzles), and 
r e (7 nozzles). 

Although one -dimensional, infinite reac- 
tion rate performance predictions can exceed 
1200 seconds at low pressure and low thrust 
levels, kinetic effects on I sp can decrement the 
equilibrium value by up to 15%. 

Kinetic losses are minimized in the high p c 
regime. Two-dimensional and boundary’ losses 
combine to produce a relatively constant 2% 1^ 
decrement for all p c and F. Maximum I sp values 
range between 930 seconds at T c of 2700 K to 
1030 seconds at T c of 3100 K. 

Although they deliver the same amount of 
thrust, multiple nozzle configurations under- 
perform their single nozzle counterpart by 
between 2 to 4% due to a decrease in area ratio 
and kinetic efficiency of each nozzle. However, 
the length of each nozzle in the cluster is 
shorter than the single nozzle length. A 65% 
nozzle length reductions can be achieved by a 5 
nozzle cluster as compared to a single nozzle 
delivering the same thrust level. 

An order of magnitude variation in the rates 
of the recombination reactions can cause 
performance decrements or increases of less 
than 2%. 


Concluding Remarks 

Geometric parameters which describe the 
near-throat region of the convergent-divergent 
nozzle were selected to be representative of a 
chemical rocket engine. The effect of these 
parameters on NTR nozzle performance was 
not investigated, but the results of this study 
indicate that hydrogen recombination is largely 
limited to the near throat region. Sensitivity 
studies of all these parameters would help to 
maximize the recombination of H to H 2 , maxi- 
mizing the energy release and delivered perfor- 
mance of an NTR nozzle. 

Other assumptions which were made to 
describe the nozzle geometry may have an 
effect on 1^ if they were changed. For example, 
the parabolic nozzle contour was not optimized 
according to known methods. Although nozzle 
performance was shown to be insensitive to 
area ratio at high chamber pressures, extrapola- 
tion of these results for a nozzle with exit radius 
of 1 meter instead of 5 meters might not 
provide accurate results. Nozzle lengths, 
assumed to be 80% of a 15' cone, may be 
substantially longer in reality for the high area 
ratio regimes. Finally, approximations were 
used to calculate the geometric variations due 
to multiple nozzle clusters. Differences may 
result from using realistic values to account for 
spacing between adjacent nozzles, etc. 

Boundary layer assumptions made were 
felt to be realistic and could be approximated in 
an actual nozzle. However, the sensitivity of 
nozzle performance to wall temperature needs 
to be investigated. 

Results of this study which may have 
impact on previous or existing studies or 
programs include: 

• no performance benefit can be 
expected from low chamber pres- 
sure for the design constraints 
used in this study (refer to Figure 

9), 

• results matched presented NER- 
VA derivative reactor (NDR) per- 
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formance data, where the operat- 
ing conditions were T c =2700 K, 
p c =1000 psia, F=75,000 lb f , area 
ratio = 500: 1 had a predicted per- 
formance of 925 seconds, com- 
pared to this study’s T c =2700 K, 
p c =1000psia, F= 100,000 lb f , area 
ratio = 463:1 predicted perfor- 
mance of 915 seconds, and 

• the maximum attainable specific 
impulse for a chamber tempera- 
ture of 3100 K is approximately 
1030 seconds, for T c of 2900 K, 
maximum performance is approx- 
imately 970 seconds, and forT c of 
2700 K, maximum performance is 
approximately 930 seconds. 
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APPENDIX 

FIGURES OF PARAMETRIC DATA 


The figures on the following pages graphically present the results of the computational study 
described in the report. Figures are grouped into sets of six plots: four specific impulse versus 
chamber pressure plots (one for each thrust level), one area ratio versus chamber pressure plot, and 
one nozzle length versus chamber pressure plot. Each of these six plots are repeated for each 
chamber temperature (three) and each multiple nozzle configuration (four). All curves are for the 
single nozzle configuration unless otherwise noted. 

Fig Description Page Fig Description Page 


A-l T c =2700 K F=10Klb f 1 Nozzle 16 

A-2 T c =2700 K F=30Klb f 1 Nozzle 16 

A- 3 T c =2700 K F=100Klb f 1 Nozzle 16 

A -4 T c =2700 K F=300Klb f 1 Nozzle 16 

A-5 T c =2700 K Area Ratio 1 Nozzle 16 

A-6 T c =2700 K Nozzle Length 1 Nozzle ..16 

A-7 T c =2900 K F=10Klb f 1 Nozzle 17 

A-8 T c =2900 K F=30Klb f 1 Nozzle 17 

A-9 T c =2900 K F=100Klb f 1 Nozzle 17 

A- 10 T c =2900 K F=300Klb f 1 Nozzle 17 

A-ll T c =2900 K Area Ratio 1 Nozzle 17 

A-12 T c =2900 K Nozzle Length 1 Nozzle ..17 

A-13 T C =3100 K F=10Klb f 1 Nozzle 18 

A- 14 T C =3100 K F=30Klb f 1 Nozzle 18 

A-15 T c =3100 K F=100Klb f 1 Nozzle 18 

A- 16 T C =3100 K F=300Klb f 1 Nozzle 18 

A- 17 T C =3100 K Area Ratio 1 Nozzle 18 

A-18 T C =3100 K Nozzle Length 1 Nozzle ..18 

A- 19 T c =2700 K F=10Klb f 3 Nozzles 19 

A-20 T c =2700 K F=30Klb f 3 Nozzles 19 

A-21 T c =2700 K F=100Klb f 3 Nozzles 19 

A-22 T c =2700 K F=300Klb f 3 Nozzles 19 

A-23 T c =2700 K Area Ratio 3 Nozzles 19 

A-24 T c =2700 K Nozzle Length 3 Nozzles. 19 

A-25 T c =2900 K F=10Klb f 3 Nozzles 20 

A-26 T c =2900 K F=30Klb f 3 Nozzles 20 

A-27 T c =2900 K F=100Klb f 3 Nozzles 20 

A-28 T c =2900 K F=300Klb f 3 Nozzles 20 

A-29 T c =2900 K Area Ratio 3 Nozzles 20 

A-30 T c =2900 K Nozzle Length 3 Nozzles. 20 

A-31 T C =3100K F=10Klb f 3 Nozzles 21 

A-32 T C =3100 K F=30Klb f 3 Nozzles 21 

A-33 T C =3100 K F=100Klb f 3 Nozzles 21 

A-34 T c =3 1 00 K F=300Klb f 3 Nozzles 21 

A-35 T C =3100 K Area Ratio 3 Nozzles 21 

A-36 T C =3100 K Nozzle Length 3 Nozzles. 21 


A-37 T c =2700KF=10Klb f 5 Nozzles 22 

A-38 T c =2700 K F=30Klb f 5 Nozzles 22 

A-39 T c =2700 K F=100Klb f 5 Nozzles 22 

A -40 T c =2700 K F=300Klb f 5 Nozzles 22 

A-41 T c =2700 K Area Ratio 5 Nozzles 22 

A-42 T c =2700 K Nozzle Length 5 Nozzles 22 

A -43 T c =2900 K F=10Klb f 5 Nozzles 23 

A-44 T c =2900 K F=30Klb f 5 Nozzles 23 

A -45 T c =2900KF=100Klb f 5 Nozzles 23 

A-46 T c =2900 K F=300Klb f 5 Nozzles 23 

A-47 T c =2900 K Area Ratio 5 Nozzles 23 

A-48 T c =2900 K Nozzle Length 5 Nozzles 23 

A-49 T c =3100 K F=10Klb f 5 Nozzles 24 

A-50 T c =3100KF=30Klb f 5 Nozzles 24 

A-51 T c =3100KF=100Klb f 5 Nozzles 24 

A-52 T c =3100KF=300Klb f 5 Nozzles 24 

A-53 T c =3100 K Area Ratio 5 Nozzles 24 

A-54 T c =3100 K Nozzle Length 5 Nozzles 24 

A-55 T c =2700KF=10Klb f 7 Nozzles 25 

A-56 T c =2700 K F=30Klb f 7 Nozzles 25 

A-57 T c =2700 K F=100Klb f 7 Nozzles 25 

A-58 T c =2700 K F=300Klb f 7 Nozzles 25 

A-59 T c =2700 K Area Ratio 7 Nozzles 25 

A-60 T c =2700 K Nozzle Length 7 Nozzles 25 

A-61 T c =2900 K F=10Klb f 7 Nozzles 26 

A-62 T c =2900 K F=30Klb f 7 Nozzles 26 

A-63 T c =2900KF=100Klb f 7 Nozzles 26 

A-64 T c =2900 K F=300Klb f 7 Nozzles 26 

A-65 T c =2900 K Area Ratio 7 Nozzles 26 

A-66 T c =2900 K Nozzle Length 7 Nozzles 26 

A-67 T c =3100 K F=10Klb f 7 Nozzles 27 

A-68 T c =3100 K F=30Klb f 7 Nozzles 27 

A-69 T c =3100KF=100Klb f 7 Nozzles 27 

A-70 T c =3100 K F=300Klb f 7 Nozzles 27 

A-71 T c =3 1 00 K Area Ratio 7 Nozzles 27 

A-72 T c =3100 K Nozzle Length 7 Nozzles 27 
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Aren Ratio Specific Impulse (sec) Specific Impulse (sec) 


Appendix: Figures of Parametric Data 


l sp for 10Klb f Thrust Level 



Figure A-l 


l sp for 30Klb, Thrust Level 



I for lOOKIb, Thrust Level 

sp I 



l sc for 300Klb, Thrust Level 



Chamber Pressure (psia) 

Figure A-4 


Nozzle Area Ratio 



Chamber Pressure (psia) 

Figure A-5 


Nozzle Length 



1C 10 

Chamber Pressure (psia) 

Figure A-6 
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Appendix: Figures of Parametric Data 


l sp for 1 0KIb, Thrust Level \ v for 30Klb ( Thrust Level 




l sp forlOOKIb, Thmst Level 



I for 300Klb, Thrust Level 

sp 1 


0 

0 

w 

'w' 

to 

3 

a 

E 

o 

0 

® 

a 

w 



Chamber Pressure (psia) 

Figure A-10 


Nozzle Area Ratio 


Nozzle Length 



!F=iOK3b. 
iF=3°KJb. I 
|F=IOOKJbJ 
IF=500KIp, 


Chamber Pressure (psia) 

Figure A-ll 



Chamber Pressure (psia) 

Figure A-12 
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Appendix: Figures of Parametric Data 



Figure A>13 


I for 30Klb t Level 

sp f 



\ v for lOOKIb, Level 


\ 9 for 300Klb, Level 




Nozzle Area Ratio 



Figure A-17 


lF=10Kt. | 
F=30Ktr, j 
F=100Ktb, 
MOOKk, 


Nozzle Length 



F=10K3b. | 
F=30Kl b, 1 
F=100K 3b r 
F=30oidb; 
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Area Ratio Specific Impulse (sec) Specific Impulse (sec) 


Appendix: Figures of Parametric Data 


l $p for 1 0KIb, Thrust Level l sp for 30Klb t Thrust Level 




l sp for lOOKIb, Thrust Level 


l sp for 300Klb, Thrust Level 





(. •' 

10 10 * 10 10 


Nozzle Length 


IF=10K£, , 
F=30Kir, j 
F-IOOK lbJ 
i F=300KJb, 



F=10KHb, 

F=30Kfc, 

F=100Kfc r 

F=300KIb, 


Chamber Pressure (psia) 

Figure A-23 


Chamber Pressure (psia) 

Figure A-24 
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Nozzle Length (inches) Specific Impulse (sec) Specific Impulse (sec) 


Appendix: Figures of Parametric Data 


l sp for 1 0Klb f Thrust Level l sp for 30Klb, Thmst Level 




l sp for 1 00Klb f Thrust Level 



Figure A-27 


l„ for 300Klb, Thrust Level 



Nozzle Area Ratio 



3 Nozzles 2 


IK : : , 

It It It It 

Chamber Pressure (psia) 

Figure A-29 


Nozzle Length 


e 

| 

a 

a 

< 


Chamber Pressure (psia) 

Figure A-30 
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Appendix: Figures of Parametric Data 


I for 1 0KIb, Thrust Level I for 30Klb, Level 

sp I sp l 




l s; for lOOKIb, Level 


1120 



Figure A-33 


I for 300Klb. Level 

sp 1 



c 

0 


< 


Nozzle Area Ratio 


Nozzle Length 



Chamber Pressure (psia) 

Figure A-35 


Chamber Pressure (psia) 

Figure A-36 
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Appendix: Figures of Parametric Data 


I. for 1 0KIb, Thrust Level l„ for 30Klb, Thrust Level 




I for lOOKIb, Thrust Level 

sp I 



Figure A-39 


L for 300Klb, Thrust Level 

sp t 



Nozzle Area Ratio 



F=10KJb, ' 
F=30KJb, i 

F=100 K3bJ 

F=300KIb f ' 


Nozzle Length 



Chamber Pressure (psia) 

Figure A-42 
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Appendix: Figures of Parametric Data 


l sp for 1 0KIb, Thrust Level l sp for 30Klb, Thmst Level 




I for lOOKIb, Thrust Level 

sp T 



Figure A-45 


L for 300Klb, Thrust Level 

sp l 




10 * ! 

10 1C 1C 10 


iwonb. 

jF =30KJ b, 

felOOH cy 

iF=300KJb; 


Chamber Pressure (psia) 

Figure A -47 


Nozzle Length 



Chamber Pressure (psia) 

Figure A-48 
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Appendix: Figures of Parametric Data 


I for 1 0KIb, Thrust Level I for 30Klb, Level 

sp 1 sp t 




l„ for lOOKIb, Level 



I. for 300Klb, Level 



Nozzle Area Ratio 



F=IOKt, 

F=30Kfc, 

F=100K3b, 


F=300KJb, 


Chamber Pressure (psia) 

Figure A-53 


Nozzle Length 



Chamber Pressure (psia) 

Figure A-54 
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Appendix: Figures of Parametric Data 


l sp for 1 0KIb, Thrust Level 



l sp for 30KI5, Thrust Level 



I for lOOKIb, Thrust Level 

sp I 


L for 300Klb, Thrust Level 

sp f 




Nozzle Area Ratio 



Chamber Pressure (psia) 

Figure A-59 



FslOKJb, 

F=30KJb, 

F=100idb, 

F=300KIb, 


Chamber Pressure (psia) 

Figure A-60 
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Appendix: Figures of Parametric Data 


l sp for 1 0KIb, Thrust Level l sp for 30Klb f Thrust Level 




l sp for lOOKIb, Thrust Level 



L for 300Klb f Thrust Level 

sp f 



Nozzle Area Ratio 


Nozzle Length 




Chamber Pressure (psia) 

Figure A-66 
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Appendix: Figures of Parametric Data 


l sp for lOKIb, Thrust Level 



I for 30Klb f Level 

sp 



l„ for lOOKIb, Level 



Chamber Pressure (psia) 

Figure A-69 


I for 300Klb, Level 

sp 1 



Nozzle Area Ratio 


Nozzle Length 




Chamber Pressure (psia) 

Figure A-72 
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